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MERGER HISTORIES OF GALAXY HALOS AND IMPLICATIONS FOR DISK SURVIVAL 
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ABSTRACT 

We study the merger histories of galaxy dark matter halos using a high resohition ACDM A^-body simulation. Our 
merger trees follow ^ 17,000 halos with masses Mg = (10^"'^ — lO^^)/i~^M0 at z = and track accretion events 
involving objects as small as m 2± IO^^/i^^Mq. We find that mass assembly is remarkably self-similar in m/MQ^ and 
dominated by mergers that are ~ 10% of the final halo mass. While very large mergers, m > 0.4 Mq, are quite rare, 
sizeable accretion events, m ~ 0.1 Mg, are common. Over the last ~ 10 Gyr, an overwhelming majority 95%) 
of Milky Way-sized halos with Mq = lGi^'^h~^MQ have accreted at least one object with greater total mass than the 
Milky Way disk (m > 5 x 10^°/i~^Mq), and approximately 70% have accreted an object with more than twice that 
mass (to > 10^^/i~^Mq). Our results raise serious concerns about the survival of thin-disk dominated galaxies within 
the current paradigm for galaxy formation in a ACDM universe. In order to achieve a ~ 70% disk-dominated fraction 
in Milky Way-sized ACDM halos, mergers involving m ^2x lO^^/i~^M0 objects must not destroy disks. Considering 
that most thick disks and bulges contain old stellar populations, the situation is even more restrictive: these mergers 
must not heat disks or drive gas into their centers to create young bulges. 

Subject headings: cosmology: theory — dark matter — galaxies: formation — galaxies: halos — 
methods: A^-body simulations 



1. INTRODUCTION 

In the cold dark matter (CDM) model of structure for- 
mation, dark matter hal os form via the continuous accre- 
tion of s maller syste ms (iPeebles )1982': 'Blum cnthal et al] 
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iCole et all (20081). 



of the type predicted can help explain many properties 
of the observed universe. Major mergers are believed to 
play an important r o le in sh aping th e Hub bl e sequence 
(iToomre fc Toomrd 119721 : IB arne p Il988t iHernauistI 



19931: iNaab fc BurkertI l2003l: lKhochfa.r fc Burkert 
2009: ICox et al.l l2006at iRobertson et aLri2006aib,ti 



Mailer et all 120061 : iJesseit et aP l2007t iBournaud et al.l 
2007() and tri ggering star form at ion and AGN ac- 
tivitv (iMihos fc Hern quistI Il996l : iKolatt et all Il999l: 

120061 iBarton et all 
explain the origin 
eTaLl fl996l: 



Cox et all I2006bl: IWoods et al 



20071 ) . Minor mergers may help expiam 
of thick disks (iQuinn et all 11993: l\yalker . . _ , 
Abadi et al.' '2 0031: iBrook et al I l2004t iD alcanton et all 
2005: KazantzidisleLalJ l2007l: iHavashi fc Chiba 2006) 
cause anti-truncation ([] 
produce ex tended diffuse li ght a 

galax i es (I Johnston et all Il996l: 

1999t [&mock et all l20o¥ iBuUock fc JohnstonI [20051 : 



1200711 . and 
light com ponents around 
iHelmi fc Whita 



Purcell et al.l 120071: iBell et all 120071 )" However, there is 
lingering concern that mergers are too common in CDM 
cosmologies for th in disk-dor ninated systems to survive 
(iToth fc Ostrikej[r9 92: Wvsi lMl fKormendv fc Fished 
l2005HKautsch et a!.i,2006. ). In this paper we present the 
merger statistics necessary for addressing this issue. 
The formation of disk galaxies within hydrody- 
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namic simulations in hierarchic al CDM co smologies 
has proven problematic (e.g. iNavarro fc S teinmet^ 
|2000( ). While there have been some successes m 
forming g alaxies with disks i n cosmological sim- 
ulations (lAbadi et all l2003l : |S ommer-Larsen et al.l 

20031: iBrook et al.l l2004l: IRobertson et all l2004l: 



Kaufmann et all l2007at iGovernato et al.l l2007h . the 

general problem is far from resolved. The resultant disks 
are often fairly thick and accompanied by large bulges, 
and the systems that form disks tend to have special 
merger histories. The resultant thick disk and bulge 
stars also tend towards a broad range of stellar ages, 
instead of being dominated by predominantly old stars. 
Moreover, the successes depend strongly on effective 
models that describe physics on scales far below the 
simulation resolution, which are poorly understood. 
Given the current difhculties in understanding ah initio 
disk formation, one can consider a less ambitious, but 
more well-posed question. Even if disk galaxies can 
form within CDM halos, can they survive the predicted 
merger histories? 

Unfortunately, the prevalence of thin-disk or even 
disk-dominated galaxies in the universe is difhcult 
to quantify with large observational samples. Some 
promising approaches use asymmetry vs. concentra- 
tion to define morphological type (e.g., Illbert et"an 
l2006l) . and some use a combination color and con- 
centration indicators (e.g. IChoi et all l2007t iPark et all 
|2007( ). Despite the wide range of definitions, the gen- 
eral consensus in the literature is that ~ 70% of ^ 
1 ^ ^ ~ ^ M(7) halos host disk - dominated, late-type galax- 
ies ( e .g. IWeinmaiin et al] l2006t Ivan den Bosch et al 



200^ Illbert et all 120061 : iChoi et all 120071 : iPark et al 
20071 ). We adopt this number for the sake of discus 



sion in what follows, but none of our primary results on 
merger statistics depend on this number. 

Also relevant to the discussion of galaxy merger his- 
tories is the prevalence of pure disk galaxies in the uni- 
verse, i^autsch et_al. (2006) have compiled a statistically 
meaningful sample of edge-on disk galaxies and found 
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that ~ 16% of these objects are "simple disks" with no 
observable bulge component. In principle, this statis- 
tic places strong constraints on the merger histories of 
galaxies. Moreover, a large fraction of disk galaxies with 
bulges contain pseudo-bulges, which may be the prod- 
ucts of secular processes and not the remnants of an 
early merger even t (e.g. iKormendv fc KennicuttI l2004l : 
iCaroUo et al1[200l . These cases provide further motiva- 
tion to quantify the predicted merger histories of galaxy 
halos in the favored CDM cosmology. 

Here we use a large dissipationless cosmological ACDM 
N-body simulation to track the merger histories of an 
ensemble of ^ 17, 000 halos with 2 = masses Mq = 
10^^ — 10^^/i~^Mq. We focus on halos of fixed mass at 
z = 0, and concentrate specifically on Milky Way-sized 
systems, Mq ~ IO^'^H-^Mq. We cate gorizc the accretion 
of objects as small as m ~ 10^^h~^MQ and focus on the 
infall statistics into main progenitors of z = halos as a 
function of lookback time. As discussed below, the main 
progenitor is defined to be the most massive progenitor 
of a z = halo tracked continuously back in time. 

A merger is defined here to occur when an infalling halo 
first crosses within the virial radius of the main progeni- 
tor. In most cases we do not track subhalo evolution after 
accretion. We have chosen to track mergers in this way in 
order to provide a robust prediction. An understanding 
of an accreted halo's subsequent orbital evolution and im- 
pact with the central disk region is essential for any com- 
plete understanding of galaxy merger statistics. However, 
this evolution will be sensitive to the baryonic distribu- 
tion within both the main progenitor halo and the satel- 
lites themselves. The halo merger rate we present is a 
relatively clean measure that can be used as a starting 
point for more detailed investigations of galaxy-galaxy 
encounters. Still, it is worth pointing out that for most 
of the mergers we consider, impacts with the central disk 
region should occur relatively shortly after accretion. As 
we show in the Appendix, events with m > O.IMq typ- 
ically happen at a redshift z when the main progenitor 
mass, Mzi is significantly smaller than Mq, such that 
the merger ratio is fairly large m/Mz > 0.2. Therefore, 
even ignoring the enhanced orbital decay that will be 
caused by a central disk potential, the dynamical friction 
decay times are expected to be short for these events, 
with central impacts occurring wit hin r < 3 Gyr for 
typical orbital parameters {B ovlan -Kolchin et "all l2007l : 
IZentner et all l2005t IZentner fc BuUockl l2003f l. As dis- 
cussed in conjunction with Figure 5 in §3, destruction 
times of ~ 3 Gyr are consistent with our measurements 
of subhalo evolution. 

The outline of this paper is as follows. In [J2]we discuss 
the numerical simulations used and the method of merger 
tree construction. In [J3]we present our principle results, 
which characterize the accretion mass functions of halos 
and the fraction of halos with mergers as a function of 
lookback time. In i]4]we discuss these results in reference 
to the problem of disk survival in a hierarchical universe, 
and we summarize our main conclusions in ^ 

2. SIMULATIONS 

Our simulation consists of 512'^ particles, each with 
mass nip ~ 3.16 x IO^H'^Mq, evolved within a comoving 
cubic volume of 80h~^ Mpc on a side usi ng the Adaptive 
Refinement Tree (ART) A^-body code (jKravtsov et al.l 



[T997l[200l . Thecosmolo gy is a fiat ACDM model, with 
parameters JIa/ = 1 — Q,a = 0.3, h = 0.7, and erg = 0.9. 
The simulation root computational grid consists of 512^ 
cells, which are adaptivcly refined to a maximum of eight 
levels, resulting in a peak spatial resolution of 1.2h^^ 
kpc, in comoving units. This simulation and the methods 
we use to construct merger trees have been discuss ed else- 
where ([AHioodiE^ [ISS IMchiEiEin llO^ Here 
we give a brief overview and refer the reader to those 
papers for a more complete discussion. 

Field dark matter halos and subhalos are identified us- 
ing a variant of th e bound density maxima algorithm 
(jKlvpin et al.lll999l ). A subhalo is defined as a dark mat- 
ter halo whose center is positioned within the virial ra- 
dius of another, more massive halo. Conversely, we de- 
fine a field halo to be a dark matter halo that does not 
lie within the virial radius of a larger halo. The virial 
radius is defined as the radius of a collapsed self gravi- 
tating dark matter halo within which the average density 
is Ayir times the mean density of the universe. For the 
family of flat cosmologies {ft m + = 1) the value of 
A„,> can be approximated by (jBrvan fc Normanlll998[ ): 

18^^ + 82(»^(z) - 1) - 39(a„(z) - 1)^ 

^vir — n / \ • \ ) 

Masses, AI, are defined for field halos as the mass 
enclosed within the virial radius, so that M — 
{4:TT/3)R^flmPc^vir- With thcsc definitions, the virial 
radius for halos of mass M at z = is given by: 

/ M \ 

Note that this mass definition based on a fixed overden- 
sity is largely conventional, and is traditionally used as 
a rough approximation for the radius within which the 
halos a re virialized. We re fer the reader to the recent 
work of lCuesta et al.l l|2007[ ) for further discussion of this 
issue. 

Halo masses become more difficult to define in crowded 
environments. For example, if two halos are located 
within two virial radii of each other, mass double- 
counting can become a problem. Also, subhalos can be- 
come tidally stripped if they are accreted into a larger 
halo. While the stripped material typically remains 
bound to the larger host halo, it is no longer bound to 
the smaller subhalo and should not be included in the 
subhalo's mass. In these cases, the standard virial over- 
density definitions arc not appropriate. In order to over- 
come this ambiguity, we always define a halo's radius and 
mass as the minimum of the virial mass and a "trunca- 
tion mass" - defined as the mass within the radius where 
the log-slope of the halo density profile becomes fiattcr 
than —0.5. This definition of truncation mass is a rela- 
tively standar d pratice w hen dealing with simulati ons of 
this kind (e.g. | Klvpin fc H oltzman 1997: Kravtso v et al.l 
l2004HZentner et al.ll2005f ). and we follow this convention 
to remain consistent with other work in this field. In 
practice, our field halos have masses and radii defined by 
the standard virial relations 98% of all non- subhalos). 
It is fair then to interpret our merger rates as infall rates 
into halo virial radii. The masses of objects just prior to 
infall are more likely affected by this definition, but the 
overall effect on our results is not large. As a test, we 
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Fig. 1. — Sample merger trees, for halos with Mq ~ 10^^'^/i~^Mq . Time progresses downward, with the rcdshift z printed on the left 
hand side. The bold, vertical line at the center corresponds to the main progenitor, with filled circles proportional to the radius of each 
halo. The minimum mass halo shown in this diagram has m = 10^'^h~^MQ. Solid (black) and dashed (red) lines and circles correspond to 
isolated field halos, or subhalos, respectively. The dashed (red) lines that do not merge with main progenitor represent surviving subhalos 
at ^; = 0. Left: A "typical" merger history, with a merger of mass m ~ O.lMg ~ O.bM^ at z = 0.51. Right: A halo that experiences an 
unusually large merger m ~ 0.4A/o ~ l.OMz a,t z = 0.65. 



have redone our main analysis using an (extrapolated) 
virial mass for inf ailing halos. The results on fractional 
merger rates change only at the ~ 5% level. 

In the event that a halo experiences a close pass with 
another halo — entering within the virial radius for a short 
time, then exiting the virial radius never to return — the 
two halos are considered isolated, even though one may 
lie within the virial radius of the other. Conversely, if the 
smaller halo falls back within the virial radius and the 
two halos subsequently merge together, then we continue 
to consider the smaller halo a subhalo even during the 
time when it lies outside the virial radius of its host. 
This has been referred to as the "stitching" method — 
as opposed to "snipping," which would count the abov e 
example as two separate mergers (iFakhouri fc Mal[200l . 

By constructing mass functions, we find that these halo 
catalogs are complete to a minimum mass lO^°/i~^M0. 
This allows us to measure the accretion of objects 10 
times smaller than 10^^h~^MQ halos, or objects up to 
1000 times smaller than IO^^H'^Mq halos. In all cases, 
we denote the mass of an accreted object as m. Over- 
all we have a total of 17,241 halos in our sample with 



Mo > lO"/i"^Af0, and (6642, 2479, 911, and 298) 
halos in logarithmically-spaced mass bins centered on 
logMo — (11-5, 12.0, 12.5, and 13.0) respectively, in units 
of h-^MQ. 

Our merger tree const ruction mirrors that described 
in iKravtsov et al.l ()2004f ) and uses 48 stored timesteps 
that are approximately equally spaced in expansion fac- 
tor between the current epoch a = [1 -\- z)~^ = 1.0 and 
a = 0.0443. We use standard terminologies for 'progeni- 
tors and descendant. Any halo at any timestep may have 
any number of progenitors, but a halo may have only one 
descendant — defined to be the single halo in the next 
timestep that contains the majority of this halo's mass. 
We use the terms "merger" and "accretion" interchange- 
ably to designate the infall of a smaller halo into the 
virial radius of a larger one. The term 771am progenitor is 
used to reference the most massive progenitor of a z = 
halo tracked continuously back in time. 

Throughout most of this work we present results in 
terms of absolute mass thresholds on the infalling mass 
m. Our principle statistics are quantified using the in- 
falling mass thresholds in terms of the final z = mass 
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of the main progenitor halo (e.g. m > O.IA/q). Indeed, 
many of our results are approximately self-similar with 
respect to halo mass when the infalling mass cut is de- 
fined in this scaled manner. By definition, the maxi- 
mum mass that a merging halo can have is m = O.SMq. 
Note that it is common in the literature to study the 
merger ratio of an infaUing object, m/Mz, where is 
the main progenitor mass at the redshift z, just prior 
to the merger. Here, Mz does not incorporate the mass 
m itself and therefore m/Mz has a maximum value of 
1.0. A parallel discussion that uses m/Mz is presented 
in the Appendix, but the absolute mass thresholds are 
used as our primary means to quantify merger statistics 
in the main part of this paper. We make this choice for 
two reasons. First, it is relatively easy to understand 
completeness effects using a fixed threshold in m, while 
completeness in m/Mz will vary as a function of time and 
will change from halo to halo depending on its particular 
mass accretion history. Second, an event with m/Mz ^ 1 
does not necessarily imply that the infalling object m is 
large compared to the final halo mass Mq. In order to 
be conservative, we would like to restrict ourselves to 
mergers that are large in an absolute sense compared to 
typical galaxy masses today. 

Figure [T] shows two pictorial examples of merger trees 
for halos with approximately equal z = masses Mq ~ 
1Q^^-^/i~^Mq. Time runs from top to bottom and the 
corresponding redshift for each timestep is shown to the 
left of each tree. The radii of the circles are proportional 
to the halo radius R ~ M^/^, while the lines show the 
descendent-progenitor relationship. The color and type 
of the connecting lines indicate whether the progenitor 
halo is a field halo (solid black) or a subhalo (dashed 
red). The most massive progenitor at each timestep — 
the main progenitor — is plotted in bold down the mid- 
dle. The ordering of progenitor halos in the horizontal 
direction is arbitrary. Once a halo falls within the radius 
of another halo, it becomes a subhalo and its line-type 
changes from black solid to red dashed. When subhalo 
lines connect to a black line this corresponds to a cen- 
tral subhalo merger or to a case when the subhalo has 
been stripped to the point where it is no longer identified. 
When field halos connect directly to a progenitor with- 
out becoming subhalos in the tree diagram it means that 
the subhalo is stripped or merged within the timestep 
resolution of the simulation. Halos that arc identified as 
subhalos of the main halo at z = are represented by 
the dashed-red lines that reach the bottom of the dia- 
gram without connecting to the main progenitor line. 

Note that the extent to which we can track a halo 
after it has become a subhalo, and the point at which a 
subhalo is considered "destroyed" is dependent both on 
spacing of our output epochs and mass resolution of the 
simulation. This is another reason why we count mergers 
when a halo falls within the virial radius (when the lines 
in Figure 1 change from solid-black to dashed-red) and 
not when a subhalo experiences a central merger with its 
host. 

The left diagram ( "halo 810" ) in Figure 1 shows a fairly 
typical merger history, with a merger of mass m ^ 0.1 Mq 
at z ~ 0.51. The merger ratio at the time of the merger 
was m/Mz — 0.5. The right diagram ("halo 853") shows 
a very rare type of merger history with a massive event 
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Fig. 2. — Average mass accretion histories for three bins in halo 
mass, as a function of lookback time. Each bin gives the average 
for bins of size logMo = 0.5, centered on the stated value. More 
massive halos accreted a larger fraction of their mass at late times. 

m ~ 0.4Mo at z ~ 0.65. This was a nearly equal-mass 
accretion event at the time of the merger, m/Mz — 1-0. 
Note that neither of these large mergers survive for long 
as resolved subhalos — they quickly lose mass and merge 
with the central halo. Each of these halos has two ~ 
lQ^'^h~^MQ subhalos that survive at z = 0. 



3. RESULTS 



3.1. Accretion Histories and Mass Functions 

The literature is rich with work on the cumulative mass 
accretion histories of halos as a function of redshift (e.g 
Wechsler et"alll2002t [Zhao et al.|[200l ITasitsiomi et al.l 



20Cil l Li et al.ll2007l and references therein). We begin 
by re-examining this topic for the sake of completeness. 
Figure [2] shows average main progenitor mass accretion 
histories, Mz = M{z), for halos of three characteristic 
final masses, Mq = M{z = 0). We confirm previous 
results that halo mass accretion histories are character- 
ized by an initial rapid accretion phase followed by a 
slower accretion phase, and that more massive halos ex- 
perien ce the rapid accretion phase later than less massive 
halos (jWechsler et al.ll2002l ). Milky Way-sized halos with 
Mq = lQ^^h~^ Mq will, on average, accrete half of their 
mass by z ~ 1.3, corresponding to a lookback time of 
- 8.6 Gyr. 

While Figure [2] provides some insight into when mass 
is accreted into halos, we are also interested in char- 
acterizing how this mass is accreted. (Ultimately, we 
will present merger statistics for a joint distribution of 
both time and mass ratio.) Now we investigate the 
mass function n{m) of objects larger than m that have 
merged into the main progenitor over its history. The 
solid line in Figure [3] shows n(m) averaged over ha- 
los in the Mq = IO^^/i-^Mq bin, plotted func- 
tion of m/Mp. On average. Milky Way-sized halos with 
Mq ~ 10^^h~^MQ experience ~ 1 merger with objects 
larger than m ~ 10^^h~^AlQ, and ~ 7 mergers with ob- 
jects larger than m ^ lO^'^h"^ Mq over the course of their 
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Fig. 3. — Mass functions of accreted material, with respect to the 
final halo mass. Lines show the cumulative number of mergers that 
a halo experiences with objects larger than m/Mo, integrated over 
the main progenitor's formation history. The (black) squares show 
the average for I0^^h~^ Mq halos; (red) crosses show the average 
for lQ^^h~^ Mq halos. Lines through the data points show the fits 
given by Equation |3] The upper/lower dashed lines indicate the 
~ 25%/20% of halos in the lO^^h'^MQ sample that have experi- 
enced exactly two/zero m > O.IMq merger events. Approximately 
45% of halos have exactly one m > O.IMq merger event; these sys- 
tems have mass accretion functions that resemble very closely the 
average. 

lives. 

For sonic purposes, an analytic characterization of the 
accreted mass function will be useful. We have investi- 
gated the average n{m) function for halos in the mass 
range Mq = 10"-^ - IO^^/i-^Mq and find that the shape 
of this function is remarkably similar over this range 
(smaller Mq were neglected in order to achieve a rea- 
sonable range in m/Mo). Specifically, we find that n{m) 
is well-characterized by a simple function of a; = m/Mo: 

n{> x) ^ Ax-" {0.5- xf, (3) 

with /3 = 2.3, a = 0.61, and x < 0.5 by construction. 
Interestingly, we find that the overall normalization in- 
creases monotonically with halo mass, and that the trend 
can be approximated as A{Mo) ~ 0.47 log]^Q(Mo) — 3.2, 
where Mq is in units of h~^MQ. This mass-dependent 
normalization, together with Equation [31 reproduces our 
measured n(m) functions quite well — to better than 5% 
at all m in smaller halos (Mq = 10"-^ - IO^^/i-^Mq), 
and (somewhat worse) to 15% at higher masses (Mq = 
1Q12.5 _ io13/i-1Mq). 5 

We would also like to understand the scatter in the 
accreted mass function from halo to halo at fixed Mq. 
It is not appropriate to simply describe the variation in 
n(m) at a fixed mass, because the total mass accreted 
is constrained to integrate to less than Mq. This means 
that the number of small objects accreted may be anti- 
correlated with the number of large objects accreted. 

^ The quoted errors are restricted to n > 0.05. 
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Fig. 4. — The fraction of halo's final mass A/q accreted from 
objects of mass m/Mo, integrated over time. The differential (up- 
per) and cumulative (lower) mass fractions are shown. The thick 
(black) lines show simulation results and the thin (red) lines show 
the prediction from standard Extended Press Schechter (EPS). Ac- 
cretions with m ~ (0.03 — 0.3)Mo dominate the mass buildup in 
halos of all Mq. The distributions are approximately self-similar 
with host halo mass. This is in reasonably good agreement with 
the EPS expectation, although the EPS fractions systematically 
sit above the simulation results. 

With this in mind, we provide an illustration of the scat- 
ter with the two dashed lines in Figure [31 The upper 
dashed line shows the average n{m) for the ~ 25% of 
halos that have experienced exactly two accretion events 
larger than m = O.IMq. The lower dashed line shows 
the average n{m) for the ^ 20% of halos that have ex- 
perienced exactly zero m > O.IMq accretion events. Ap- 
proximately ^ 45% of halos have exactly one m > O.IMq 
event, and these have an average accreted mass function 
that is very similar to the overall average shown by solid 
(black) line in Figure [H Halos with fewer large mergers 
show a slight tendency to have more small mergers, but 
the effect is not large. 

Figure m presents some of the same information shown 
in Figure [21 but now in terms of the mass fraction, f{m), 
accreted in objects larger than m for Mq = 10^^ and 
10^^h~^MQ halos (thick lines, see legend). The up- 
per panel in Figure [3l shows the differential fraction, 
d//dlnm = (— TO^/Afo)dn/dm, while the lower panel 
plots the integrated fraction f{m). As before, we have 
normalized the accreted masses, to, by the final 2 = 
main progenitor mass Mq. We find that /(> m) is 
also well fit by Equation [3] (to better than 10% across 
aU masses Mq = lO^^ _ IO^'^H-^Mq)'^. As before, 
a; = m/Mo, but now A{Mq) ~ 0.17 logio(Mo) - 0.36, 
with Mq still in units of H-^Mq. The best fit param- 
eters are a = 0.05, and /3 = 2.3. The lines are trun- 
cated at m/Mo = 0.01 and 0.001, corresponding to our 
fixed resolution limit at to = 10^^h~^ Mq. The thin 
(red) lines of the same line types show the same quan- 

^ The quoted errors are restricted to / > 0.1. 
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Fig. 5. — Left: The fraction of Milky Way-sized halos, Mq ~ 10 h ^Mq, that have experienced at least one merger larger than a 
given mass threshold, m, since look-back time t. Right: The three solid (black) lines show the fraction of Mq ~ 10^^ Mq halos that 
have experienced at least one, two, or three mergers larger than m = O.IMq in a lookback time t. The upper solid line is the same as the 
solid lino in the left panel. The dashed (red) line include only mergers that are "destroyed" before 2 = (i.e. lose > 80% of their mass 
by z = 0). Conversely, the dot-dashed (blue) line shows the fraction of halos with at least one merger that "survives" as a subhalo as a 
function of lookback time to the surviving halo's merger. Large accretions of m > O.IMq typically do not survive for more than ~ 3 Gyr 
after accretion. 



titles predicted fro m Extended Press Schechter (EPS 
iLacev fc Cold[l993h Monte-Carlo merger trees. Each of 
these lines is based on 50 00 trees generated using the 
iSomerville fc KolattI (|1999D algorithm. 

In broad terms, the mass spectrum of accreted ob- 
jects agrees fairly well with the EPS expectations, espe- 
cially considering the relative ambiguity associated with 
defin i ng halo masses i n simulations (e.g. ICohn fc White! 
[20071: iDiemand etall [20071 : iCuesta et all [200l . How- 
ever, it is worth discussing the similarities and dif- 
ferences in some detail. It is a well-known expec- 
tation from EPS that the total mass accreted into 
a halo of mas s Mp is dominated by objects of mass 
m ~ G.lMn (ILacev fc Cold 119931 [Zentner fc Bullock! 
[200l!PurceU et al.!!2007! : !Zentneii!2007! ). Our simulations 
reveal that indeed m ~ (0.03 — 0.3)Mo objects are the 
most important contributors to the final halo mass. 

EPS trees predict self-similar mass fractions across all 
halo masses. Our more massive halos, however, show 
a slight tendency to have more of their mass accreted 
in collapsed objects, across all scaled masses m/Afg- As 
discussed in association with Equation[3]above, the over- 
all normalization of the mass spectrum is slightly higher 
for our more massive halos. For example, the mass ac- 
creted in objects larger than m = 0.01i\jfo is ^ 45% for 
Mo = 10^2 halos and ~ 50% for Mq = 10" halos. Both 
of these fractions are low compared to the ~ 65% ex- 
pected from the semi-analytic EPS merger trees. 

It is interesting to estimate the total mass fraction ac- 
creted in collapsed objects (m > 0) by extrapolating the 
n{m) fit given in Equation [3| to m — > 0. Wc find 



^0.5 

Ax^^" (0.5-a;)'^dx 
/o 

^ 0.24 A 



(4) 



/(>0) = 



0.5A/o 



dm 



dm 



In the second step, P and a are the shape parameters 
in Equation [3] In the last step we have used our best 
fit parameters a = 0.61 and P = 2.3. Using Equation 2 
for A{Mo), we find that the total accreted mass fraction 
(in virialized halos) increases from / ~ 0.50 to 0.70 as 
Mo varies from 10"-^ to lO^/i-^Mg. This suggests that 
a significant fraction (~ 30 — 50%) of dark halo mass is 
accreted in the form of "diffuse" , unvirialized material, 
and that smaller halos have a higher fraction of their 
mass accreted in this diffuse form. Of course, our conclu- 
sion on "diffuse" accretion may be due, at least in part, 
to difficulties in precisely defining halo masses — and, in 
particular, in defining halo masses in dense environments 
where mergers are occurring rapidly. It is also possible 
that the accreted mass function steepens below our res- 
olution limit, resulting in a lower diffuse fraction than 
we expect from our extrapolation, but there is no clear 
physical reason to expect such a steepening. 

3.2. Merger Statistics 

Understanding how galaxy mergers can affect galaxy 
transformations and morphological fractions necessarily 
requires an understanding of halo merger statistics. Of 
specific interest is the overall fraction of halos that have 
had mergers within a given look back time. 

The left panel of Figure [5| shows the fraction of Milky 
Way-sized halos (Mq 10^^/i~^Mq) that have experi- 
enced at least one "large" merger within the last t Gyr. 
The different line types correspond to different absolute 
mass cuts on the accreted halo, from m > 0.05 A/o to 
m > 0.4 A/q. The tendency for lines to flatten at large 
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Fig. 6. — The fraction of halos of mass Mq at z = that have experienced a merger with an object more massive than O.IMq (left) and 
O.SMq (right) in the last t Gyr. Error bars are Poissonian based on the number of halos used in each mass bin. Note the fairly weak mass 
dependence. 



lookback times is a physical effect and is not an arti- 
fact of limited resolution. Specifically, the lines flatten 
at high z because the halo main progenitor masses, Af^, 
become smaller than the mass threshold on m (see Fig- 
ure 2). We find that while less than ~ 10% of Milky 
Way-sized halos have ever experienced a merger with 
an object large enough to host a sizeable disk galaxy, 
(m > 0.4 Afo ~ AxIO^^H^^Mq), an overwhelming major- 
ity ('^ 95%) have accreted an object more massive than 
the Milky Way's disk (m > O.OSMq IO^^H-^Mq). 
Approximately 70% of halos have accreted an object 
larger than m ^ 10^-^h~-^ Mq in the last 10 Gyr. 

While the left panel of Figure [5] illustrates the frac- 
tion of halos with at least one merger in a given time, 
the right panel of Figure [5] provides statistics for multi- 
ple mergers. We again focus on Mq lO^^h"^ Mq ha- 
los, but restrict our statistics to accretions with m > 
0.1 Mo ~ 10"/i^^A/o. The upper solid (black) line 
shows the fraction of halos with at least one accretion 
event within the last t Gyr. This reproduces the solid 
(black) line in the left-hand panel, but now the ver- 
tical scale is logarithmic. The middle and lower solid 
(black) lines in the right panel show the fraction of halos 
with at least two and at least three mergers larger than 
0.1 Afo in the past t Gyr, respectively. Roughly ~ 30% 
of Milky Way-sized halos have experienced at least two 
accretion events larger than ~ 10^^h~^MQ. Multiple 
events of thi s kind could be i mportant in forming ellipti- 
cal galaxies ('HernqTiist"1993':'Bovlan-Kol chin et al.ll2()05l : 
[Robertson et al. 2006a; Naab ct al. 2006). 

As discussed above, our merger events are defined at 
the time when the accreted halos first cross within the 
virial radius of the main progenitor. This definition al- 
lows us to focus on a robust statistic that is less likely to 
be affected by baryonic components. In SJH we speculate 
on the implications of these results for disk stability. In 
this context, one may be concerned that some fraction 



of our identified "mergers" will never interact with the 
central disk region, but instead remain bound as surviv- 
ing halo substructure. We expect this effect to be most 
important for smaller accretions, as larger merger events 
will decay very quickly. 

The dashed (red) line in Figure [5] shows a statistic that 
is analogous to the upper black line — the fraction of ha- 
los that have had at least one merger in the last t Gyr 
— except we have now restricted the analysis to include 
only objects that are "destroyed" before z = 0. We define 
an object to be "destroyed" if it loses more than 80% of 
the mass it had at accretion because of interactions with 
the central halo potential. (Our results are unchanged 
if we use 70 — 90% thresholds for mass loss). Likewise, 
the dot-dashed (blue) line shows the same statistic re- 
stricted to "surviving" objects. Only ^ 10% of Milky 
Way-sized halos have surviving massive substructures at 
z ~ that are remnants of m ~ 10^^h~^MQ accretion 
events. These survivors were typically accreted within 
the last ~ 3 Gyr. 

We conclude that large accretions that happened more 
than ~ 3 Gyr ago would have had significant interac- 
tions with the central disk regions of the main progen- 
itor. Indeed, this must be the case if the Milky Way 
is "typical" . The most massive Milky Way satellite, 
the LMG, was likely accret ed with a mass no larger 
than m ^ 4 X 10^°/i ~^M(d (Ivan der Marel et all [200^ : 
[Robertson et al.|[2005D . According to Figure[31 we expect 
that the Milky Way has accreted at least ^ 3 objects that 
are larger than the LMC over its history. The expecta- 
tion is that most of the larger objects that were accreted 
have been shredded by the central galaxy potential, and 
have deposited their stars in the extended stellar halo 
(jBuUock fc Johnstonl |2005( ) . Interestingly, a recent re- 
analysis of the LMG's motion suggests that it is indeed 
on it s first passage about the Milky Way (jBesla et al.l 
[2007I) . as we would expect for surviving, massive satel- 
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lites. 

Given that galaxy morphology fractions are observed 
to change syste matically as a fu nction of luminosity or 
mass scale (e.g. iPark et al.ir2007l ). we are also interested 
in exploring merger fractions as a function of Mq. The 
left panel of Figure [6] shows the fraction of halos that 
have had a merger larger than m — 0.1 Mq within the 
last t = 2, 4, ... 12 Gyr, as indicated by the t labels. The 
right panel shows the same statistic computed for larger 
m > 0.3 Mq accretion events. The most striking result 
is that the merger fraction is fairly independent of Mq. 
This suggests that halo merger statistics alone cannot 
explain the tendency for early-type spheroidal galaxies 
to reside in more massive halos. Baryon physics must 
instead play the primary role in setting this trend. (The 
right-panel in Figure [6] does show evidence that halos 
larger than Mq ^ 3 x 10^^h~^MQ have a higher fraction 
of recent {t < 2 Gyr) large (m > 0.3Mo) events, but the 
overall merger fraction is small (< 10%), even for group- 
sized halos with Mg ~ 10^^h~^MQ.) When counting 
instantaneous merger rates there appears to be a rela- 
tively universal merger r ate across a wide rang e of host 
masses Mq - lO^^-lO" (jFakhouri fc Mai [200% . Given 
this uniformity, we speculate that the weak trend we see 
in Figure [6] may ultimately be a direct result of the trend 
for more massive halos to form later, as shown in Figure 

m 

If a merger occurs early enough, we might expect 
the main progenitor halo to grow significantly after the 
merger. The fraction of mass accreted since the last 
merger is potentially important, as, for example, it could 
enable the regrowth of a destroyed disk. One could rea- 
son that the fraction of halos that experience a large 
merger and then subsequently fail to accrete a signifi- 
cant amount of mass is the most relevant statistic for 
evaluating the probability of disk formation. However, 
we find that in most cases, very little mass is accreted 
after an event that is large relative to Mq. 

Figure [7] shows the average fraction (AM/Mq) of a 
halo's final mass Mq that is accreted since the last large 
merger. Each curve corresponds to a different thresh- 
old in m. As in Figure [6l the overall trend with final 
mass Mq is very weak. We define the accreted mass by 
AM = MQ—Mm, where Mm is the mass of the host halo's 
main progenitor after the most recent large merger. The 
upper (solid, black) line includes all halos that have had 
at least one event larger than m = 0.1 Mq and the low- 
est line (long-dashed, red) includes all halos that have 
had at least one merger larger than m = 0.35 Mq. This 
shows that the fractional mass accreted since a merger 
is a decreasing function of m/Mo. It also shows that the 
fraction is typically smaU, with AM/Mq < 30% (10%) 
for m > 0.1 Mq (0.3 Mq) mergers. These trends are a 
consequence of the fact that we are doing this calcula- 
tion at a fixed Mq — if a merger m is large compared 
to Mq, then there is little room in the mass budget for 
new material to be accreted after the merger. This im- 
plies, for example, that if a disk is destroyed as a result 
of a large-m/Mo accretion event, it is unlikely that a 
new "disk-dominated" system can be regrown from ma- 
terial that is accreted into the host halo after the merger. 
However, gaseous ma terial involved in the merger may re- 
form a disk (see e.g. IZurek et al.lll988t [Robertson et al.l 
l2006aD . 



4. DISCUSSION 
4.1. Milky Way Comparison 

The Milky Way has a dark ma tter halo of mass Mq ~ 
IO^^/i-^Mq (jKlvpin et al.l 120021 ) . and its stellar mass 
is dominated by a thin disk of mass ^ 3.5 x IO-'^'^Mq 
(jKlvpin et al1[200l : lwidrow fc Dubinskil[2005h . The thin 
disk has vertical sc ale height that is j ust ^ 10% of its 
radial scale length (ISiegel et all l2002l : lJuric et all l2005l : 
iNewbe rg et al."2006V and contains stars as old as ~ 10 
Gyr (jNordstrom ct al. 2004). Moreover, stars in the lo- 
cal thick disk arc predominantly older than ~ 10 Gyr, 
and the bulge is old as well. This suggests that there 
was not significant merger activity in the Milky Way to 
drive gas towar ds the bulge or to thicken the disk in the 
past ~ 10 Gyr (lWvsdl2001h . 

Based on our results, a galaxy like the Milky Way 
would seem rare in a ACDM universe. Roughly 70% 
of dark matter halos of mass Afo — lO^'^h~^MQ have ex- 
perienced a merger with a halo of mass IO^^H'^Mq in 
the past ~ 10 Gyr. A merger of this size should thicken 
the existing disk and d rive gas into the center of the 
galaxy to create a bulge (jKazantzidis et al]|2007f) . If the 
Milky Way has not experienced a merger of this mag- 
nitude, that would make our galaxy a rare occurrence 
(< 30% of halos). On the other hand, if the Milky Way 
has experienced such a merger, it is difficult to under- 
stand its observed early-type morphology and thin-disk 
properties. 

4.2. Morphological Fractions and Thick Disks 

The degree to which the Milky Way halo is typical for 
its mass is becoming better understood thanks to the 
advent of large, homogeneous astronomical sky surveys. 
As mentioned in the introduction, broad-brush catego- 
rizations of "late type" vs. "early type" suggest that 
~ 70 % of Milky Way-sized halos host late-type galaxies 
/e.g. IWeinmann et all l2006t Ivan den Bosch et all l2007l : 



lllbert et al.l 120061 iChoi et all 120071 : iPark et all I2OO7I) . 

The degree to which "late-type" is synonymous with 
"thin disk-dominated" is difficult to quantify with cur- 
rent data sets, but for the sake of this discussion, we 
will assume that this is the case. Also discussed earlier 
were the results of Kautsch et al. (2006), who found that 
~ 16% of disk galaxies are bulgeless systems. This sug- 
gests that ~ (0.7)(0.16) - 11% of Milky Way-sized halos 
host pure disk galaxies. 

The observed morphological fractions may be com- 
pared to the halo merger fractions presented in Figure 
O These results show that an overwhelming majority of 
Milky Way-sized halos (~ 95%) experience at least one 
merger larger than the current mass of the Milky Way 
disk (> 5 X 10^°/i"^Mq). Figure |3| shows that a typical 
Mq ~ lO^^/i^^M0 halo has merged with ~ 2 — 3 objects 
of this size over its history. It is possible that merg- 
ers of this characteristic mass are responsi ble for creat- 
ing thick disk components in most g alaxies () Walker et al.l 
119961 iDalcanton fc Bernsteinll2002t i. More detailed sim- 
ulations will be required to test whether disks are de- 
stroyed or overly thickened by the predicted infall of 
m ~ 5 X lO^°/i~^M0 objects, and whether these thicken- 
ing events happen too late to explain thick dis ks as old 
as those observed (jPalcanton fc Bernstein|[2002l ) . Under- 
standing how bulgeless galaxies could exist in halos with 
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Fig. 7. — The average fractional change in a halo's main progen- 
itor mass AM /Mo since its last large merger, shown as a function 
of the halo's z = mass Mq. Each line includes only halos that 
have had a merger larger than the m value indicated {i.e. only 
halos that have had a m > 0.35Mo accretion event are included in 
the lowest line). Error bars shown are Poisson on the number of 
halos used in the average. 

mergers of this kind is a more difRcult puzzle. 

Perhaps more disturbing for the survival of thin disks 
are the statistics of more substantial merger events. 
Figure shows that m > 10^^h~^MQ accretions are 
quite common in Milky Way-sized halos, with ~ 70% 
of Mq ~ 10^^h~^ AIq objects experiencing such a merger 
in the past ~ 10 Gyr. Of course, the impact that these 
events will have on a central disk will depend on orbital 
properties, gas fractions, and star formation in the merg- 
ing systems. Generally, however, a merger with an object 
~ 4 times as massive as the Milky Way thin disk would 
seem problematic for its survival. 

We find that a small fraction of a halo's final mass is 
typically accreted into the main progenitor subsequent 
to m ~ O.lMo mergers — this suggests that the re- 
growth of a dominant disk from material accreted af- 
ter such a merger will be difficult (see Figure [71). We 
conclude that if ~ 70% of Milky Way-sized halos con- 
tain disk-dominated galaxies, and if the adopted ACDM 
cosmology is the correct one, then mergers involving 
m 2± IO^^H^^Mq objects must not result in the destruc- 
tion of galaxy disks. This is a fairly conservative con- 
clusion because if we naively match the percentages of 
mergers with an early- type fraction of ~ 30%, then Fig- 
ure 5 suggest that the critical mass scale for disk survival 
is significantly larger, to > 2 x 10^^h~^AlQ. Specifically, 
mergers involving objects that are 5 times the current 
Milky Way disk mass must not (always) destroy disks. 

We remind the reader that the lookback times depicted 
in Figures 5 and 6 correspond specifically to the times 
when infalling halos first fall within the main progeni- 
tor's virial radius. Our estimates suggest that the corre- 
sponding central impacts should occur ^3 Gyr later for 
the mass ratios we consider. Therefore, when we quote 



merger fractions to a lookback time of ~10 Gyr, this will 
correspond to an actual impact ~7 Gyr ago. Of course, 
the infalling systems will also lose mass as they fall to- 
wards the central galaxy. As we have emphasized, the 
detailed evolution of merging objects can only be deter- 
mined with focused simulations, and the outcome of the 
subsequent mergers will depend on the baryonic compo- 
nents and orbital properties of the systems involved. 

Kazantzidis et al. (2007) have performed a focused 
N-body simulation in order to investigate the morpho- 
logical response of a thin, Milky Way type stellar disk 
galaxy to a series of impacts with 6 satellite halos ^ of 
mass TO ~ (l-2)xlOi°M0 (~ 30-60% of the disk mass). 
They find that a dominant "thin" stellar disk component 
survives the bombardment, although its scale- height was 
seen to increase from 250 pc to ~ 400 pc, and a sec- 
ond ~ 1.5 kpc "thick" component was also created. In 
addition, a new central bar / bulge component was also 
generated in these fairly small encounters. While it is 
encouraging to see that a thin disk can survive some 
bombardment, mergers with objects ~ 6 times as mas- 
sive as those considered by Kazantzidis et al. should be 
very common in Milky Way-sized halos. It remains to be 
seen how the infall of m ^ IO^^H'^Mq objects will af- 
fect the morphologies of thin, ^ 4 x IO^^Mq stellar disks. 
Moreover, the merger history considered by Kazantzidis 
et al. extended to fairly recent events. While there was 
no explicit star formation prescription in these simula- 
tions, one would expect a broad range of stellar ages in 
the thickened disk stars in this case, instead of a pre- 
dominantly old population — as seems to be observed in 
actual galaxies. 

4.3. Morphology-Luminosity Trends 

Another well-established observational tr end is the 
morphology-luminosit y relation (recently, IChoi et al.l 
120071: iPark et al.ll2007l) . which, when interpreted in terms 
of a morphology-halo mass relation, demonstrates that 
the fraction of late-type galaxies contained within dark 
matte r halos is anti-correlat ed with the mass of the 
halo (jWeinmann et al] I2006D . For large galaxy ha- 
los, Mq — 10^^h~^MQ, the late-type fraction is just 
~ 30 %, compared to ^ 70% for Milky Way-sized sys- 
tems (IWeininann et al.ll2006l: Ivan den Bosch et al.ll2007l : 
lllbert et al.l 120061 IChoi et all 120071 : IPark et all 120071). 
The result presented in Figure 6 is perhaps surprising 
in light of this fact. Specifically, merger histories of 
galaxy halos are almost self-similar in Mq when the in- 
falling mass m is selected to be a fixed fraction of Mq. 
For example, - 18% (70%) of Mq ~ IO^^H-^Mq halos 
have experienced an to > 0.3Mo (O.IMq) merger event 
in the last 10 Gyr. This fraction grows only marginally 
to ~ 25% (80%) for Mq - IO^^H-^Mq halos. The im- 
plication is that dark matter halo merger histories alone 
cannot explain the observed correlation between early- 
type fraction and halo mass. Specifically, baryon physics 
must play the primary role in setting the observed trend 
between galaxy morphology and halo mass. 

We note that their 6 accretions were chosen from a high- 
resolution N-body simulation, and that this number of events is 
fairly typical of what we find based on ~ 2500 Milky Way-sized 
halos in Figure 4. However, it is also typical of a Milky Way-size 
halo to experience 1 — 2 mergers more massive than any of these 
accretions. 
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4.4. Successive Minor Mergers 

Recently, Bournaud ct al. (2007) have used focused 
simulations to investigate the response of a very massive 
~ 2 X IO^^A/q disk-dominated galaxy within a lO^^M© 
halo to mergers with total mass ratios ranging from 
m/Mz ~ 0.02 to 1.0. Broadly speaking, they find that 
m/Mz =0.1 merger ratio events can transform their 
disk galaxy to an SO, and that m/Mz = 0.3 ratio events 
produce ellipticals. It is unclear how these ratios would 
change for a smaller primary disk mass, considering that 
their disk mass is extremely massive for a halo of this 
size. In comparison, ^ 95% of our Milky Way-sized halos 
experience an event with a merger ratio of m/Mz > 0.1 
(corresponding to m > 0.05 Mq, see Appendix) in the 
last 10 Gyr. Similarly, ^ 60% of our halos experience 
m/Mz > 0.3 events (See Figure IS]). We note that the 
results of this type of simulation will be sensitive to the 

:as fractions an d ISM model of the interacting galaxies 

Robertson et al. 20 06a). 

iHavashi fc Chibal (|2006f ) have also investigated the re- 
sponse of a galactic disk to a succession of minor mergers 
of CDM subhalos. They find that subhalos more mas- 
sive than 15% of the disk mass must not merge into the 
thin disk itself, or it will become thicker than the ob- 
served disk of the Milky Way. While our merger rates 
are for subhalos entering the virial radius of the halo, 
not when it penetrates the disk, we expect the fraction 
of disk mergers involving objects of mass 6 x lO^Modot 
to be quite high. Recall that 95% of halos experience an 
accretion event larger than 5 x lO^^Afodot- Even if these 
halos lose 90% of their mass before disk impact (which 
seems unlikely) they still meet the Hayashi and Chiba cri- 
terion. Note, however, that a detailed thin/thick decom- 
positi on may be required in o rder to fully evaluate this 
limit (jKazantzidis et ani2007| ). A more detailed study of 
large mergers, including the necessary baryon physics, is 
required to fully explore this issue. 

4.5. Gas-rich Mergers 

Many of our results provide qualitative support to the 
idea that cool gas-fractions play a fundamental role in 
governing the morphological outcome of la r ge mergers 
jRobertson et al.ll2006at Brook et al.]l2007al lbl- ICox et all 
120071) . with gas-rich mergers essential to the formation 
and survival of disk galaxies. While dark halo merger his- 
tories are approximately self-similar in Mq, gas fractions 
are known to decreas e systematically with halo mass (at 
least a t z = 0, see e.g. lGeha et al.|[2006t iKaufmann et al.l 
l2007bl and references therein). This implies that gas-rich 
mergers should be more common in small halos than in 
large halos. If gas-rich mergers do allow for the forma- 
tion or survival of disk galaxies, then the gas-fraction- 
mass trend may provide an important ingredient in ex- 
plaining the observed morphology-mass trend. Specifi- 
cally, small halos should experience more gas-rich merg- 
ers, while large halos should experience more gas-poor 
mergers. However, stars resulting from a gas-rich merger 
will be younger than the lookback time of the merger, 
suggesting that mergers of this type only serve as an 
adequate explanation for early mergers. Also, these gas- 
rich mergers would require sufficient angular momentum 
to keep the resulting disk from being too centrally con- 
centrated. 



4.6. Gomparison to Previous Work 
Semi-analytic mode ls (e.g. IDiaferio et al.l l2001l : 



lOkamoto fc Nagashimal |2001D and models base d on 
cosmological SPH simulations (|Maller et al.ll2006[ ) have 
demonstrated that many of the observed trends between 
galaxy morphology, density, and luminosity can be 
explained in ACDM-based models. However, these 
results rely on the ability to choose a characteristic 
galaxy mass ratio for morphological transformation 
'"char (usually Tchar 0.3). The assumption is that 
that galaxy-galaxy mergers with mgai/Afgai < 7'char 
allow a disk to remain intact, while all lar ger mergers 
produ ce spheroids. In a recent investigation. iKoda et al.l 
(|2007t l have explored a two-parameter model, where 
stellar spheroid formation depends on both m/Mz and 
the absolute halo mass, Mz- Host halos with masses 
smaller than a critical Mz are assumed to have very low 
baryon content in this picture. Using the PINOCCHIO 
density-field algorithm (.Monaco ct al...200^ to generate 
halo merger statistics, iKoda et al.l ( 20071 ) find that the 
fraction of disk-dominated galaxies can be explained if 
the only events that lead to central spheroid formation 
have m/Mz > 0.3 and Af;, > 4 x W^H'^Mq. 

The goal of this work has not been to find a ratio that 
can explain morphological fractions. Rather, our aim has 
been to emphasize the relatively large number of more 
minor mergers that could potentially be ruinous to disk 
survival. As discussed in the Appendix, events with m ~ 
O.lAfo in galaxy halos typically have m/Mz — 0.2 at the 
tmie of accretion. These ~ lO"ft."^Af0 mergers would 
produce no morphological response i n the primary disk 
under many staiid ard treatments (e.g. lMaller et al.l[2006l : 
iKoda et"aLll2007l ). For the sake of comparison. Figure [S] 
in the Appendix shows galaxy halo merger statistics for 
fixed m/Mz cuts. Our results are in good agreement with 
those derived by Koda et al. (2007; Figure 4) usi ng the 
PINO CCHIO algorithm, and with those quoted bv lWvsd 
(|2006[ ) for an analysis made by L. Hebb using the GIF 
simulations. 

Finally, we mention that ICole et all ()2007[ ) have used 
the large Millennium Simulation to investigate the pro- 
genitor mass functions of halos. In qualitative agreement 
with our findings. Cole et al. (2007) find that the fraction 
of mass coming from halo progenitors is lower than ex- 
pected from standard EPS treatments. Their approach 
was somewhat different than ours, as they focused on 
the full progenitor mass function, as opposed to the mass 
function of objects that merged into the main progenitor 
as we have here. They estimate that ^ 14% of a halo's 
progenitors are not accounted for in collapsed objects 
at any rcdshift. This may be compared to our estimate 
of ~ 30 — 50% for the fraction of mass not directly ac- 
creted in the form of virialized objects for Mq = 10^^ to 
IO^^H'^Mq halos. The differences between our numbers 
and theirs may come from the fact that we have actually 
measured slightly different quantities. We also used dif- 
ferent halo-finding algorithms and halo mass definitions, 
and utilized different formulations to extrapolate the sim- 
ulation results to unresolved masses (a peak heights for- 
mulation in their case, and a direct mass function formu- 
lation in our case). A more thorough investigation of the 
differences associated with halo-finding algorithms and 
mass definitions is reserved for a future paper. For the 
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purposes of this work, it is useful to point out that while 
a direct comparison is difhcult to make at this time, if 
anything, our results on over all merger c o unts seem low 
compared to the results given iCole et alJ (pOOl . 

5. CONCLUSION 

We have used a high-resolution ACDM iV-body 
simulation to investigate the merger histories of ~ 
17,000 galaxy dark matter halos with masses Mq ~ 
2q11-13^-i^^ at z = 0. Mergers with objects as small 
as m = 10^^h~^ Mq were tracked. The principle goal has 
been to present the raw statistics necessary for tackling 
the issue of thin disk survival in ACDM and for provid- 
ing a cosmological context for more focused simulations 
aimed at understanding the role of mergers for processes 
like morphological transformation, star formation trig- 
gering, and AGN fueling. 

Our main results may be summarized as follows: 

1. Mass accretion into halos of mass Mq at z = is 
dominated by mergers with objects of mass m ~ 
(0.03-0.3)Mo (FigureHl). Typically, - 1-4 merg- 
ers of this size occur over a halo's history (Figure 
[3]). Because these mergers tend to occur when the 
main progenitor's mass, Mz, was somewhat smaller 
than Mq, these dominant events have fairly large 
merger ratios, m/Mz — 0.1 — 0.6 (see Appendix). 

2. The mass accretion function, n{m), of mergers 
larger than m accreted over a halo's history is well- 
described, on average, by a simple analytic form, 
n{x = m/Mo) = Ax"" (0.5 - xY , with a = 0.61 
and P = 2.3. The normalization increases as a 
function of the halo's mass at z = 0, Mq, as 
A ~ 0.471ogio(Mo) - 3.2. By extrapolating this 
fit, we find that the total mass fraction accreted 
in objects of any mass (to > 0) does not asymp- 
tote to 1.0, but rather increases with AIq from 
~ 50% in Mo = 10" '5/i-iMo halos to ~ 70% in 
Mq = 10^^h~'^MQ halos. This suggests that a non- 
zero fraction of a halos mass may be accreted as 
truly "diffuse" material. 

3. An overwhelming majority (95%) of Milky Way- 
sized halos with Mq ~ 10^'^h~^ Mq have accreted 
an object larger than the Milky Way's disk (to > 
5 X lO^°ft."^M0) in the last 10 Gyr. Approximately 
70% have had accretions with to > 10^^h~^MQ 
objects over the same period, and 40% have had 
TO > 2 X lO"/i"^M0 events (Figure [^j. 

4. Halo merger histories are approximately self- 
similar in to/A/q for halos with masses in the range 
Mo = 10" - IO^^H-'^Mq (FigureE]). This suggests 
that the empirical trend for late-type galaxies to 
be more common in smaller halos is not governed 
by differences in merger histories, but rather is as- 
sociated with baryon physics. 

5. Typically, a small fraction, 20 — 30%, of a halo's 
final mass AIq is accreted after the most recent large 
merger with to > (0.1 - 0.2) Mq objects (Figure[3). 
This suggests that the "regrowth" of a disk from 
newly accreted material after a large merger is un- 
likely. Note that this does not rule out the pos- 



sibility that a disk reforms from gaseous material 
involved in the merger itself. 

The relatively high fraction of halos with large to ~ 
0.1 Afo merger events raises concerns about the survival of 
thin disk galaxies within the current paradigm for galaxy 
formation in a ACDM universe. If we naively match per- 
centages using Figure m we find that in order to achieve 
a - 70% disk-dominated fraction in Afo = lO^^/i-^Afg 
halos, then to ~ 0.2 Afo ~ 2 x lO^/i-^Afo ~ 3 x lO^Afg 
objects must not (always) destroy disks. Furthermore, 
since stars in the local thick disk and bulge are predom- 
inantly older than ~ 10 Gyr, this suggests that these 
mergers in the past ~ 10 Gyr must not drive gas towards 
the bulge or significantly thicken the disk. Note that the 
total mass in such an accreted object is ~ 10 times that 
of the Milky Way disk itself. Moreover halos typically do 
not accrete a significant fraction of their final mass after 
these mergers (~ 20% on average). Finally, as noted in 
the Appendix, to 0.2 Mo events typically have merger 
ratios of m/Mz — 0.4 at the time of the merger. These 
numbers do not seem encouraging for disk survival, and 
may point to a serious problem with our current under- 
standing of galaxy formation in a ACDM universe. 

Our basic conclusion is unlikely to be sensitive to un- 
certain cosmological parameters. Note that the simula- 
tions considered here have a fairly high erg = 0.9. At a 
fixed f2m , a lower cts will syste matically produce slightl y 
more recent merger events (e.g. lZentner fc Bullockl2003( ). 
However, given that the merger fractions measured using 
m/Mo are approximately self-similar in Mo (and there- 
fore in Mo/Af^), we expect that the overall merger frac- 
tions will be fairly insensitive to power-spectrum normal- 
ization. 

As discussed in the introduction, a complete investi- 
gation into the issue of disk survival will require an un- 
derstanding of the orbital evolution of objects once they 
have fallen within the main progenitor halo's virial ra- 
dius and on the subsequent impact of interacting galax- 
ies. Both of these outcomes will depend sensitively on the 
baryonic components in the main halo and in the smaller 
merging object. For this reason, the present work has fo- 
cused on halo mergers, defined to occur when an infalling 
halo first crosses within the virial radius of the main pro- 
genitor halo. The merger statistics presented here are 
relatively devoid of uncertainties and can be used as a 
starting point for direct simulations of galaxy-galaxy en- 
counters. Simulations of this kind will be essential to 
fully address the broader implications of these frequent, 
large mergers, which seem to pose a serious challenge to 
disk survival. 

The simulation used in this paper was run on the 
Columbia machine at NASA Ames. We would like to 
thank Anatoly Klypin for running the simulation and 
making it available to us. We are also indebted to Bran- 
don AUgood for providing the merger trees. We thank 
Chris Brook, T. J. Cox, Fabio Governato, Christian 
Maulbetsch, Brant Robertson, Matias Steinmetz, Jerry 
Sellwood, and Rosemary Wyse for helpful comments on 
an earlier draft, and Onsi Fakhouri and Chung- Pci Ma for 
useful discussions about halo merger rates. We thank the 
anonymous referee for several insightful suggestions that 
helped to improve the quality and clarity of the paper. 
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APPENDIX 

The bulk of the paper focused on statistics of halo mergers at a fixed absolute mass threshold m on the merging 
objects. In this section we present statistics for m/Mz — the merger ratio relative to the main progenitor mass Mz 
at the rcdshift z prior to accretion. By definition, Mz < Mq, causing merger statistics to show larger accretions 
when presented in terms of m/Mz- For example, the merger ratio equivalent to Figure [3] (bottom panel) fits to 
/(> x) ~ Ax~"(l - x)f^ instead of Equation[3] (to better than 10% across mass range 10^^-^ - lO^^'^/i^^M©), with 
parameters A{Mq) ~ 0.05 logxQ(A/o) -2.0, a = 0.04, and f3 = 2.0, and Mq in units ofh~^MQ. The m/Mz equivalent for 
n{> x) also fits to this equation (to better than 10% in the mass range 10^^ ° — Ky^^-^h'^^Mo) with parameters a ~ 0.4, 
(3 ~ 2.0, and A{Mq) ~ 2.0 log]^Q(il/o) — 21.5. Given that our simulation output has a fixed physical mass resolution, 
rrires = lQ^^h~^MQ, and that halo main progenitor masses, Mz, vary as a function redshift, the completeness in m/Mz 
will, in principle, vary from halo to halo as a function of each particular mass accretion track and lookback time. We 
have accounted for this variable completeness limit in an average sense in what follows. 

The left panel of Figure [8] shows the fraction of Mq = 10^^/i~^Mo halos that have experienced at least one merger 
larger than a given threshold ratio {rt ~ 0.1, 0.2, ... 0.9) in m/Mz in the last t Gyr. This result may be compared 
to Figure [5l which shows the analogous fraction computed using fixed absolute m cuts relative to Afo- In order to 
provide results that are robust to our completeness limit in m, for each threshold cut, rt, the lines are truncated at 
the lookback time when the average Mq ~ lQ^'^h~^ Mq halo's mass falls below Mz = TOics/?'*- We see that ^ 50% 




Fig. 8. — Left: The fraction of Milky Way-sized halos, Mq ~ 10^'^h~^ Mq, that have experienced at least one merger larger than a given 
ratio m/Mz since look-back time t, where Mz is the main progenitor mass at the time of accretion. Lines are truncated at epochs where 
the mass resolution of the simulation limits our ability to resolve a given m/M^ ratio (i.e. when Mz gets so small that a quoted m/Mz 
ratio falls below the resolution with m > 10^^h~^ Mq). Right: The average change in mass AM/Mq that a halo experiences since its last 
major merger. Different lines show different "major merger" ratios, m/Mz, where the ratio is defined relative to the main progenitor mass 
at the time of accretion. Error bars show Poissonian errors on the number of host halos averaged. 
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Fig. 9. — Le/t; The fraction of halos of mass Mq at 2 = that have experienced a merger ratio larger than m/Mz = 0.3 in the last 
t Gyr. The ratio is defined relative to the main progenitor mass at the time prior to accretion Mz- Error bars are Poissonian based on 
the number of halos used in each mass bin. Note the strong mass trend, in contrast to the results presented in Figure |6] where mergers 
were defined on a strict mass threshold, rather than merger ratio. Right: The fraction of halos of mass Mq at z = that have experienced 
a merger ratio larger than m/Mz = 0.6 in the last t Gyr. The ratio is defined relative to the main progenitor mass at the time prior to 
accretion Mz. Error bars are Poissonian based on the number of halos used in each mass bin. Note the strong mass trend, in contrast to 
the results presented in Figure |6] where mergers were defined on a strict mass threshold, rather than merger ratio. 

of halos have had a m/Mz > 0.4 event in the last ~ 10 Gyr, and that ^ 10% have experienced nearl y eq ual-mass 
mergers, m/Mz > 0.9, in this time. These results are in good agreement with similar results quoted bv IWv sc (200^ 
for an analysis made by L. Hebb using the GIF simulations. 

The right panel of Figure [S] shows that, typically, the largest m/Mz events occur before most of the final halo mass 
Mq is accreted. Specifically we show the fraction of mass AM/Mq accreted since the last major merger. Results for 
different major merger ratio thresholds are shown as different line types. Each line presents the average AM/Mq at 
fixed Mq and we only include halos that have actually had a merger of a given ratio (within the last ^11 Gyr) in this 
figure. Among halos that have had nearly equal-mass merger events, m/Mz ^ 0.9, the fraction of mass that is accreted 
since that time is significant, with AM/Mq ~ 70% for Mq = IO^^H'^Mq halos. If we associate post-merger accretion 
with the potential "regrowth" of a galactic disk, then, by comparison with Figure [3 high merger-ratio events are less 
of a concern for disk formation than high m/Mp ratio events. Unlike most high m/Mz events, mergers that arc large 
relative to Mq typically have very little (< 20%) fractional mass accretion after the merger. 

The left panel of Figure [5] shows the fraction of halos that have had a merger larger than m/Mz = 0.3 within the 
last i = 2, 4, ... 10 Gyr, as a function of halo mass Mq. The right panel shows the same statistic computed for larger 
m/Mz > 0.6 mergers. Unlike the result shown in Figure [6l there is a fairly significant mass trend, with more massive 
halos more likely to have experienced a major merger at a fixed lookback time. Note, however, that most of these 
"major mergers" involve relatively small objects in an absolute sense. Much of this trend is driven by the fact that 
Mz falls off more rapidly with z for high AIq halos compared to low Mq halos. 

The previous discussions lead to explore the relationship between an accreted halo's absolute mass m and the mass- 
ratio it had when it was accreted, m/Mz- The two panels of Figure [TUl illustrate this relationship for objects accreted 
into Mq = lO^^/i~^M0 halos. The thick, solid line in the left panel shows the median merger ratio, m/Mz^ that a 
halo of mass to/Mq has when it was accreted. Specifically, we plot m/Mz given m/MQ in this diagram. We see that 
typically m/Mz — 2m/Mo, as shown by the dashed line in the figure. The dotted lines show the 68 % spread in the 
distribution of m/Mz given m/Mo- The opposite relationship is shown in the left panel of Figure [TOl Here, the thick, 
solid line shows the median m/Mg value given a merger of mass ratio m/Mz- We see that the majority of high-ratio 
events involve objects that are small compared to the final halo mass Mq. Note that this result, and the associated 
distributions, are complete only to mergers that occur within the past ~ 11 Gyr. If we were able to track main 
progenitor masses Mz back to arbitrarily early times, we would expect a very larger number of high m/Mz events with 
small absolute m/Mo values. 
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Fig. 10. — Left: An illustration of P(m/Mz\m, Mo) — the distribution of merger ratio m/Mz at the time of accretion given a value of 
m for Mq = IO^^/i^^Mq halos. We plot m/Mz vs. m/Mo for clarity. The solid line shows the median and the dotted lines show the 68% 
spread. Right: An illustration of P{m/Mo\m/Mz) — the distribution of merging masses m given a merger ratio m/Mz- We see that the 
majority of high m/Mz events occur with Mz is small compared to the final halo mass Mq. Therefore most high-mass ratio mergers are 
small m mergers in an absolute sense. 



